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WATER-LEAVING  RADIANCE  OF  DEEP  AND  SHALLOW 
WIND-ROUGHENED  WATER  BODY  ILLUMINATED  BY  SUN 
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Naval  Research  Laboratory,  Ocean  Optics  Section,  Code  7333, 
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e-mail:  <vihaltrin@nrlssc.navy.mil,  Telephone:  228-688-4528 

1.  INTRODUCTION 

This  paper  is  devoted  to  derive  remote  sensing  coefficient  (RSC)  or  radiance  reflection  coefficient  is  a 
principal  product  of  atmospheric  correction  algorithms  applied  to  the  remotely  measured  optical  images  of  the  ocean. 
The  RSC  contains  information  about  angular  structure  of  light  radiance,  roughness  of  the  ocean  surface,  and  optical 
properties  of  the  water.  This  presentation  analyses  RSC  and  presents  it  as  a  product  of  three  physically  different 
values  or  coefficients.  The  first  coefficient  depends  on  geometrical  parameters  of  illumination  and  atmospheric 
optical  properties,  the  second  coefficient  depends  on  roughness  of  the  ocean  surface,  and  the  third  coefficient 
depends  on  the  optical  properties  of  seawater.  The  approach  used  to  derive  the  third  in-water  coefficient  is  valid  for 
all  levels  of  water  turbidity  -  from  the  clearest  open  waters  to  the  most  turbid  waters  of  river  mouths  and  extremely 
scattering  waters  of  Yellow  Sea  type.  It  is  shown  that  presented  structure  of  RSC  allows  to  increase  precision  of 
retrieving  seawater  optical  properties  from  remotely  measured  optical  images  of  the  ocean. 

2.  REMOTE  SENSING  COEFFICIENT 

By  definition,  remote  sensing  coefficient  (RSC)  of  the  sea  or  ocean,  or  brightness/radiance  coefficient  of  the 
sea  or  ocean  (as  it  is  defined  in  European  literature),  is  a  ratio  of  upward  radiance  Lu  to  downward  irradiance  Ed  of 
light  just  above  the  water  surface, 

r„  =  L'M-  (') 


Equation  (1)  could  be  rewritten  as 


L+  L+  E~  U  TtDT, 

p(6)=  -M- 

H  C'd  nw 


Pie), 


(2) 


where  Lu  is  upward  radiance  just  below  the  water  surface,  Ed  is  upward  irradiance  just  below  the  sea  surface, 
T°  =  Euj Lu  —  7]°  Ini  is  a  transmittance  of  diffuse  light  ascending  from  the  water  by  the  wavy  surface,  7j°  is  a 
transmittance  of  atmospheric  diffuse  light  by  the  wavy  surface,  nw  is  the  water  index  of  refraction,  7j,  =  Ed  /  Ed  is 

a  transmittance  of  total  atmospheric  light  by  wavy  sea  surface,  and  p{6)  =  Lu  j  Ed  is  radiance  coefficient  of  water 
body  of  finite  depth,  and  0  is  a  zenith  angle.  Here  superscripts  (+)  and  (-)  denote,  respectively,  values  above  and 
below  the  water  surface. 

The  total  transmittance  of  water  surface  can  be  expressed  as  follows, 

n=/r,+o-/)[(i-9.)n°+4n1.  o) 

where  /  is  a  portion  of  the  water  surface  covered  by  foam,  Tj  =  1  —  Aj  is  a  foam  transmittance,  Aj  is  a  foam 

albedo,  qa  is  a  portion  of  downward  irradiance  in  atmosphere  just  above  the  water  surface  that  is  due  to  a  direct  solar 

s  s  s 

light,  Ti  =  1  —  Rf  is  a  transmittance  of  direct  solar  light  by  the  sea  surface  free  of  foam,  and  RF  is  a  Fresnel 

reflection  coefficient  of  direct  solar  light  by  wavy  water  surface  uncovered  by  foam. 

20060130  295 
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The  total  radiance  (or  brightness)  coefficient  of  water  body  p{6)  is  connected  with  the  total  diffuse  reflection 
coefficient  r=eJe~  by  the  following  relationship  [1], 


R  =  2n  jQp(d)cos6smdd9 ,  0<R<1. 


For  Lambertian  reflector,  p{6)  =  p  =  const  and  R  =  71  p. 

The  total  radiance  coefficient  p(6)  can  be  split  into  two  parts,  the  radiance  coefficient  for  diffuse  light  pD , 
and  radiance  coefficient  for  direct  solar  light  ps, 

p(,0)  =  (\-qw)pD(d)  +  qwps(6),  (5 

where  qw  is  a  portion  of  downward  irradiance  in  water  just  below  the  water  surface  [1-3].  The  value 

q,  =  e;/(ez+e;  )  can  be  expressed  through  the  similar  value  qa  as  follows.  The  downward  diffuse  and  direct 
irradiance  just  below  the  water  surface  are, 


EWD  =  f(l- Af)(EaD  +  Eas)  +  (\-  ,  Ews  =(l-f)Eas  Tf 


consequently, 


qw  =  (l-f)j-qa- 


By  combining  Eqs.  (2),  (5)  and  (7)  we  have, 

D 

rr.=it{Tip°-q.Q- ms  [p„-ft(0)]} 

nw 

td  (* 

s  -M/(1-  Af) +(i-/)[o-?,)7-°+?<  r/  t,s(9)]}k, 

here  we  assumed  that  pD  =  R/k  and  denoted  7]s(0)  =  ps(6)/ pD  .  The  value  f]s  is  dependent  on  a  zenith  angle 
6.  and  inherent  optical  properties  of  seawater  The  diffuse  transmittance  is  a  transmittance  of  direct  light  Tj5 
averaged  over  zenith  angle  6  with  the  weight  proportional  to  the  sky  radiance: 


n°=(n!L. 


so,  Eq.  (8)  can  be  rewritten  as  follows: 


■{/(l-  Af)  +  (l-/)[(l-^)(r/)s  +q„  Tf  (Js(9)]}«. 


where  R  is  a  diffuse  reflection  coefficient  of  the  water  body  illuminated  by  diffuse  light. 

According  to  [2,  3]  solar  to  sky  light  ratio  parameter  qa  is  only  a  function  of  atmospheric  parameters, 


Eg  Ta  1 


%a-  1  + 
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where  Ta  is  a  total  atmospheric  optical  thickness,  Ba  is  a  probability  of  backscattering  of  light  in  atmosphere,  and 
jls  =  COsZs ,  with  Zs  as  a  solar  zenith  angle.  The  foam  fraction  /  is  a  function  of  a  wind  speed.  According  to  [1] 
it  is  expressed  through  the  wind  speed  u  as  follows: 


f  = 


(  1.2- 10“  V’3, 

[l.210_5«3-3  (0.221m-  0.99), 


u  <9m  /sec, 
u  >  9m /sec. 


(12) 


The  transmittance  of  direct  light  by  wavy  sea  surface  is  expressed  through  the  Fresnel  reflection  coefficient  of  a 
flat  surface  and  wind  speed  as  follows  [4]: 

Tf(Zs,u)  =  1-  a0(u)~  R°(Zs){al(u)+  R°F(Zs)\a2{u)+  a3(u)Rp(Zs)^ ,  0<  w<12m/sec,  (13) 

here 

<(zs)=\ 

is  a  Fresnel  reflection  coefficient  of  directed  light  from  clear  water,  and  coefficients  at  (i  =  0,...  3)  are  given  by  the 


following  equations: 

c^{u)  =  0.001(6.944831  -1.912076  m  +0.03654833  u),  (15) 

(\{u)  =  0.743 1368  +0.0679787  u  -  0.0007 171m2,  (16) 

c^iu)  =  0.5650262  +0.0061502w- 0.0239810  w2  +0.0010695  u,  (17) 

Qg  (m)  =  -0.4128083-0.1271037m+0.0283907m2 -0.0011706m3.  (18) 


The  equation  for  transmittance  of  diffuse  light  by  wavy  sea  surface  according  to  Eqs.  (9)  and  (13)  is, 

Ti  («)  =  (Ti  )e  =  1  ~  «0 («)  -  (i?F){a1(M)  +  ( RF)[a2(u )  +  a3  (u)  (/?F)]} ,  0<  u  <12  m  /  sec,  (19) 

where 

1  rnl2 

(rf)=2  l  Rf(0)mde,  J0  F(6)de=  1,  (20) 

here  F(9)  is  a  normalized  radiance  distribution  of  the  sky  (without  including  the  direct  sunlight).  For  Lambertian 
sky  F(0)=  2/tt  we  have  the  following  equation  for  diffuse  transmittance: 

T?(u)  =  1.367 -10-5(46.434-  m)(1410  +20.6m  +u2),  0<M<12m/sec.  (2la) 

For  overcast  sky  F(6)=  2(1  +  2 COS  9)f  (4  +  7C)  (cardioid  distribution)  we  have  a  different  equation  for  diffuse 
transmittance: 

Ti(u)  =  6.123 -10-6  (59.3- u)(2564  +  33.74m  +  u),  0<«<12m/sec.  (2lb) 


cos  Zs  - , Jn2w  -  sin 2  Zs  'l  f  n2  cosZs  -  -Jnl  -  sin2Z, 


ycos  Zs  + 


4*1  - si 


sin  Z 


s  ) 


\2' 


nlcosZs  +  4n2w  -sin2Z5  J 


(14) 
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The  difference  between  values  of  transmittance  given  by  Eq.  (21a)  and  (21.b)  is  about  10%.  The  average  of 

transmittance  values  given  by  both  equations  is  about  (r°(«>).  “0.9.  The  more  elaborate  procedure  to  average 

Fresnel  reflection  coefficient  should  include  an  actual  radiance  distribution  of  the  sky  which  is  dependent  on 
atmospheric  optical  properties. 


3.  NONLINEARITY  OF  APPARENT  OPTICAL  PROPERTIES  IN  RESPECT  TO  GORDON’S 

PARAMETER 


Remote  sensing  coefficient,  as  we  show  below,  is  proportional  to  diffuse  reflection  coefficient.  Both  these 
values  are  apparent  optical  properties,  that  is  they  are  dependent  not  only  on  inherent  optical  properties,  but  also  on 
conditions  of  illumination  and  geometrical  and  optical  structure  of  the  surface  and  bottom  of  the  sea. 

Diffuse  reflection  coefficient  (DRC)  of  water  body  is  an  informative  part  of  remote  sensing  reflectance  [5]  of 
light  by  the  ocean.  DRC  contains  information  on  content  of  dissolved  and  suspended  substances  in  seawater.  DRC  is 
an  apparent  optical  property  that  depends  not  only  on  inherent  optical  properties  of  the  seawater,  but  also  on  the 
parameters  of  illumination.  The  dependence  on  inherent  optical  properties  is  expressed  through  the  dependence  on 
Gordon’s  parameter,  i.e.  the  ratio  of  backscattering  coefficient  bB  to  the  sum  of  absorption  a  and  backscattering 
coefficient,  g  —  bB  /  (a +  bB).  In  the  open  ocean  DRC  is  linearly  proportional  to  g.  This  linear  equation  is  very 
good  for  the  Type  I  open  ocean  waters  [6],  It  is  also  acceptable  for  about  90%  of  coastal  waters.  Theoretical  and 
numerical  analysis  show  that  the  linear  relationship  can  be  adequately  used  only  when  Gordon’s  parameter  g  is 

relatively  small,  i.e.  g<  0.1.  This  criterion  is  always  satisfied  in  open  ocean  waters.  The  available  database  of 
experimental  measurements 
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Figure  1.  Distribution  of  occurrences  of  measured  Gordon’s 
parameter  in  Yellow  Sea  in  2000  [7], 


Figure  2.  Dependence  of  diffuse  reflectance  coefficient 
errors  on  Gordon’s  parameteter  [7], 


show  that  in  coastal  waters  Gordon’s  parameter  may  exceed  this  critical  value  of  0.1.  In  some  very  turbid  coastal 
waters  it  can  even  reach  values  higher  than  0.95.  For  example,  in  waters  of  Yellow  Sea  or  coastal  waters  close  to 
river  estuaries  the  percentage  of  cases  when  g  >  0. 1  can  reach  50%  or  more. 

Until  recently  we  had  no  significant  cases  of  in  situ  measurements  that  show  existence  of  sea  waters  with  high 
cases  of  Gordon’s  parameter  g  exceeding  0.1.  Relatively  recent  measurements  made  in  2000  by  NRL  researches 
with  collaboration  of  Korean  scientists  in  Yellow  Sea  shows  that  these  highly  scattering  waters  exhibit  these 
characteristics  [7],  Figure  1  shows  a  histogram  of  the  frequency  of  measurements  of  Gordon’s  parameter  in  this 
expedition.  In  this  case  more  than  50%  of  measurements  correspond  to  Gordon’s  parameter  exceeding  critical  limit 
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of  linear  approach.  It  means  that  processing  of  such  data  should  involve  nonlinear  equations  connecting  DRC  with 
g.  For  that  reason  we  will  use  a  self-consistent  radiative  transfer  approach  that  takes  into  account  non-linearity  of 
inherent  optical  properties  on  Gordon’s  parameter  [8-10]. 


4.  RADIANCE  REFLECTANCE  RATIO  AND  DIFFUSE  REFLECTANCE  OF  WATER 


According  to  the  theory  of  Ref.  [8-10]  we  can  write  the  radiance  reflectance  ratio  as: 


where 


ls(0)  = 


Ps(0) 

Pd 


2(1  +  k-cos0)  [/c  —  ln(l  +  k*)]  ’ 


(22) 


K  = 


m-p2) 

1  +  Ji2 


P  = 


t 


1  -8 

+  2g+ilg(4  +  5g) 


bB 

a  +  bB  ’ 


(23) 


a  is  an  absorption  coefficient  of  water,  bB  is  a  backscattering  coefficient  of  light  by  water,  g  is  a  Gordon’s 
parameter,  and  has  physical  meaning  of  average  cosine  of  diffuse  light  inside  water,  and  6  is  a  zenith  angle 

measured  inside  water  body.  If  we  want  to  express  Eq.  (22)  through  a  zenith  angle  $  outside  water,  we  should  use 
the  Snellius  law, 


Figure  3.  Dependence  of  parameter  Tfs  on  viewing  angle 
for  different  values  of  Gordon’s  parameter. 


Figure  4.  Dependence  of  parameter  Tls  on  Gordon’s 
parameteter  for  different  viewing  angles. 


sin$  =  n^sin#,  which  results  in  substitution  in  Eq.  (22):  COS0  — >  -J\  —  sin2 1?/ .  Function  ris(0)  is 
normalized  as  follows: 

rir/2 

2jo  r]s(6)cosdsinddd  =  l.  (24) 


34 


Proceedings  of  the  III  International  Conference  “Current  Problems  in  Optics  of  Natural  Waters,”  ONW’2005,  Eds.  Iosif  Levin 

and  Gary  Gilbert,  St.  Petersburg,  Russia,  2005 


The  dependence  of  factor  T]s  on  viewing  angle  and  Gordon’s  parameter  is  shown  in  Fig.  3  and  4. 

The  diffuse  reflectance  of  shallow  sea  of  arbitrary  turbidity  and  absorption  is  represented  by  the  following 
equation  [8]: 


R^l-AbR^  +  jAb-RJe-^ 
(l-AbR0)  +  (Ab-RJR0e-v*>  ’ 


(25) 


where 


R  = 


yl+flj 


Ro  = 


r2  +J[^ 

2-Jt 


R„,  v=2  n(a  +  bB), 


(26) 


Ab  is  an  albedo  of  a  bottom,  and  zb  is  a  bottom  depth,  and  is  a  diffuse  reflectance  of  optically  deep  sea 
illuminated  by  diffuse  light. 

Equations  (22)  and  (25)  take  into  account  multiple  scattering  of  light  inside  water  and  multiple  reflection  of 
light  from  the  bottom.  It  is  valid  for  arbitrary  depths  (0  <  zb  <°°),  arbitrary  bottom  albedo  0  <  Ah  <1,  and 
arbitrary  inherent  optical  properties  of  water,  0  <  a  <  °°,  and  0  <  bB<co. 


5.  CONCLUSION 


The  result  of  this  paper  is  an  equation  that  expresses  remote  sensing  reflection  coefficient  as  a 
function  of  twelve  parameters: 

Ks(Ba  5  W,  ,Zb,  Ab,X)  — 

=  {/(“)U  -  A/A)]  +  [1  -  /(»)][{1  -  T„)}  +  (27) 

nnw 

+qa(Ba’'Ca)Ti  (Zs,u)  ris(a,bB,$ )]  }R(a,bB,zb,  Ab), 

where  Ba  is  a  light  backscattering  probability  in  atmosphere, 

Ba  =  ^  2pa(cos6)smdd6 1 \o  pa(cos9)sin0d6  (28) 

(here  pa  is  an  atmospheric  light  scattering  phase  function,  and  6  is  a  scattering  angle),  Ta  is  a  total  optical 
thickness  of  atmosphere,  u  is  a  wind  speed  at  the  level  of  water  surface,  Aj  is  an  albedo  of  a  wind  generated  foam 
or  whitecaps,  Zs  is  a  solar  zenith  angle,  tlw  is  a  refraction  index  of  water,  I?  is  a  viewing  zenith  angle  of  an  optical 
receiver,  a  is  an  absorption  coefficient  of  light  in  water,  bB  is  a  backscattering  coefficient  of  light  in  seawater,  Zb  is 
a  water  depth,  Ab  is  a  bottom  albedo,  and  X  is  a  wavelength  of  light.  Equation  (27)  shows  spectral  (wavelength) 
dependence  of  RSC  only  through  the  foam  albedo.  In  reality,  the  following  parameters  also  have  spectral  dependence 
on  the  wavelength  of  light:  absorption  and  backscattering  coefficients  of  water,  bottom  albedo,  backscattering 
probability  and  total  optical  thickness  of  atmosphere,  and,  to  some  extent,  refraction  index  of  water.  These 
dependencies  are  not  shown  in  the  right  side  of  Eq.  (27)  in  order  to  make  it  less  bulky.  The  transmittivity  of  diffuse 
light  by  the  wavy  water  surface  7j°  is  expressed  by  Eqs.  (21),  or  can  be  calculated  from  Eq.  (19)  using  current 
atmospheric  properties.  The  share  of  water  surface  covered  by  foam  /  is  computed  using  Eq.  (12).  The  share  of 
direct  sunlight  in  the  total  light  radiance  on  the  water  level  qa  is  computed  using  Eq.  (11).  The  transmittance  of 
direct  light  by  the  wavy  surface  is  computed  with  Eqs.  (13)-(18).  The  ratio  of  direct  light  reflectance  to  the 
reflectance  of  diffuse  light  T]s  is  given  by  Eqs.  (22)-(23).  And,  finally,  the  diffuse  reflectance  coefficient  of  diffuse 
light  by  the  shallow  water  body  is  given  by  Eqs.  (25)-(26).  Equation  (27)  is  derived  to  be  valid  in  the  whole  range  of 
variabilities  of  optical  properties  of  the  water  body  and  atmosphere  and  for  arbitrary  depth  of  the  water.  It  is  intended 
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to  be  used  in  algorithms  of  remote  processing  of  optical  data,  including  derivation  of  water  optical  properties,  water 

depth  and  bottom  albedo. 
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